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At strong pump powers, a semiconductor optical cavity passes through a Hopf bifurcation and 
undergoes self-oscillation. We simulate this device using semiclassical Langevin equations and assess 
the effect of quantum fluctuations on the dynamics. Below threshold, the cavity acts as a phase- 
insensitive linear amplifier, with noise ~ 5x larger than the Caves bound. Above threshold, the 
limit cycle acts as an analog memory, and the phase diffusion is ~ 10 x larger than the bound set by 
the standard quantum limit. We also simulate entrainment of this oscillator and propose an optical 
Ising machine and classical CNOT gate based on the effect. 


Many problems in simulation, optimization and ma¬ 
chine learning are analog in nature and mapping them 
onto a digital processor incurs significant overhead. As a 
result, there has been a recent revival of interest in ana¬ 
log or “neuromorphic” computing systems mn]. Devices 
that can spontaneously oscillate are a key component in 
this neuromorphic architecture. Such devices can func¬ 
tion as an analog memory [2] , a phase-insensitive ampli¬ 
fier [3l 11] , or a complex-valued neuron [5] , among other 
things. In addition, large networks of such oscillators can 
be applied to complex optimization and machine learning 
tasks, such as Ising problems [T]. 

In most dynamical systems, spontaneous oscillations 
arise from a Hopf bifurcation [6]. In optics, the sim¬ 
plest such system is the non-degenerate optical para¬ 
metric oscillator (OPO), which behaves as a quantum- 
limited amplifier below threshold [7] and has a symmet¬ 
ric limit cycle above [8|. In addition, cavity quantum 
electrodynamics (QED) systems can self-oscillate in the 
right conditions 0 E]. However, nanofabrication with 
materials such as KTP and LiNbOs is still in its 
infancy m, and most implementations of cavity QED - 
trapped atoms, quantum dots, NV centers - are not scal¬ 
able with current technology. To realize neuromorphic 
computing with photonics, there is an unfulfilled need 
for self-oscillating photonic devices based on a scalable 
technology. 

Eree-carrier dispersion can fulfill this unmet need. This 
effect is present in silicon and all HI-V semiconductors, 
and is scalable and low-power m- Previous work by 
Malaguti et al. [laiis] and Chen et al. El showed that 
when the photon and carrier lifetime are comparable, an 
optical cavity can pass through a Hopf bifurcation and 
undergo self-oscillation. However, these studies focused 
on the many-photon classical limit, where quantum fluc¬ 
tuations can be ignored. If such a device is optimized 
for low power, quantum fluctuations in the photon and 
carrier number may substantially alter the dynamics and 
limit the performance of real devices. 


* Electronic address: rhamerly@stanford.edu 


In our previous work m, we derived a set of stochastic 
equations for free-carrier optical cavities that model these 
quantum fluctuations. Here, we apply those equations 
to study the effects of quantum noise on the free-carrier 
Hopf bifurcation. 

Sections |I] and [H] discuss the general theory of the oscil¬ 
lations, which arise from an instability in the linearized 
model around the system’s fixed point. Because this is 
done in a general, scale-invariant way, it should be pos¬ 
sible to observe these oscillations in a wide range of sys¬ 
tems spanning orders of magnitude in speed, size and 
energy. Next, we consider the equations of motion close 
to the bifurcation point and show that the bifurcation re¬ 
sembles the non-degenerate OPO at threshold with some 
extra noise. Section models the device below thresh¬ 
old: it functions as a phase-insensitive linear amplifier 
with noise ^ 5x above the Caves bound m- The near¬ 
threshold behavior, which follows the critical exponents 
of the Hopf bifurcation, is discussed in Section |TV| 

The above-threshold case is covered in Section [Vl Like 
the non-degenerate OPO, the free-carrier cavity has a 
limit cycle in this regime. The above-threshold OPO 
can be considered a “quantum-optimal” limit cycle in 
the sense that it can function as an optimal homodyne 
detector. By comparison, the free-carrier limit cycle is 
^ 10X noisier than the OPO. This difference is due to 
the incoherent nature of carrier excitation and decay. 

Limit-cycle devices can be very useful in optimization 
and machine learning. In Section [VI A[ we propose and 
simulate an Ising machine based on the free-carrier limit 
cycle, which should be several orders of magnitude faster 
and less power-consuming than a supercomputer. In ad¬ 
dition, Section |VIB| discusses an all-optical XOR gate 
based on the limit-cycle effect. 


I. CONDITIONS FOR SELF-OSCILLATION 
A. Equations of Motion 

A single-mode free-carrier optical cavity has three de¬ 
grees of freedom: two field quadratures (a, a*) and the 
free carrier number N. Typically, the following effects 
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are relevant: 

1 . Cavity-waveguide coupling. This gives rise to a lin¬ 
ear loss K in the cavity field. 

2. Linear and two-photon absorption. The former 
dominates for near-bandgap operation of direct-gap 
semiconductors; the latter for indirect-gap systems. 
Gives rise to a linear loss term 77 and a quadratic 
loss term p. Both act as source terms for the carrier 
number. 

3. Free-carrier dispersion / absorption. The cavity de¬ 

tuning shifts as a function of the carrier number: 
A A - 1 - ScN. li dc = — i 62 is complex, this 

accounts for free-carrier absorption as well. 

4. Carrier decay. Typically due to recombination at 
surface sites or diffusion out of the cavity. This 
gives rise to a linear loss term 7 for N. 

In this text, we ignore the following effects: 

1 . Excitons, which tend to be the dominant effect only 
at low temperatures or in exotic materials. 

2. Thermo-optic effect. Temperature changes much 
more slowly than the photon or carrier number, so 


does not typically play a role in the fast dynamics of 
the device. It may, however, lead to stability issues, 
which are not the focus of this paper p!8l[2Q] . 

3. Optomechanical effects, which are negligible unless 
a cavity has been specifically engineered to probe 
them. 

Under these assumptions, the device can be modeled 
as an open quantum system that couples to a Markovian 
bath; see generally [2TII23] . The full quantum theory is 
quite involved and is discussed in our previous paper m- 
In short, starting from a quantum model with a bosonic 
photon mode and many fermionic carrier modes, one can 
construct a generalized Wigner function in terms of a 
set of bosonized operators and derive a Fokker-Planck 
equation for this function using the truncated Wigner 
method [SUES]- This can be recast as a set of stochas¬ 
tic differential equations (SDE’s) which sample from the 
Wigner function as a probability distribution. Assuming 
that dephasing and thermalization are much faster than 
the photon or carrier lifetimes, one obtains the following 
stochastic equations of motion (Eqs. (C16-17) from [E]): 


da = 




— (/3 -h — i(A + N 6 c) 


adt - + -y/ridf3rj - 2 ^/Pa'"dPp - \/ 2 N 52 dpfca ( 1 ) 

( 2 ) 


dU 


dN = [r]a*a + ^{a*af - dt + v^(aV/3^ +adp;) + ^/^{{a*fdpp + a^d^p)*) + ^f^Ndw^ 


d^N 


and the output optical field is: 


dPout = — dt + dPu 


(3) 


In these equations, d/din is a complex Wiener process representing the input field, which for vacuum input has 
the Ito rule d^i^dP"^^ = dt/2. The processes dp^ and dfdfca correspond to linear, two-photon and free-carrier 
absorption respectively, and also have vacuum statistics. The dw^ is a real Wiener process satisfying dre^ = dt, giving 
the Poisson statistics of carrier decay. The real and imaginary parts of 5c are 5c = 5i — 2 ^ 2 . Typical values for the 
parameters in dip are given in Table | 

These equations resemble the coupled-mode equations used to analyze semiconductor microcavities elsewhere in 
the literature [THU]- Unlike the equations used elsewhere, ( p!p| ) include quantum-noise terms. As a result, these 
equations allow us to model the quantum behavior of devices previously only discussed classically, and study the 
fundamental quantum limits to device performance. 

We can analyze optical bistability and self-oscillation by linearizing these equations of motion about their equilibrium 
point. Defining the doubled-up vector x = {5a,5a*,5N)^ the equations of motion take the following form: 


5a 


5a* 

= 

5N 



dx 


+ 


7 ^ - i(A + N5c) - 2{I3 -h ix)(a*(a 


(77 - 1 - 2/3(a*(a)(a* 


— —i5ca 

- i(A + N5c) - 2{I3 -h ix)<^*<^)* {-i5cay 
(77 - 1 - 2f3a*a)a —7 


Ax dt 


-\/k 0 

0 -a/k 

0 0 


d/din 

dPL 


+ 


-y/qdpn- 2y^a*dl3i3 - y '2NS2 dpfca 

- 2Xpadldl - VW^d/3}^^ 

y/ri{a*djdn + adp*) + y/P{{a*)‘^djdp + a^^dPi^)*) + X^dw^^ 


Bd/3i, 


5a 

5a* 

5N 


dt 


(4) 


Fdw 































3 


Likewise, the output can be related to the input and 
internal state by: 


By the cyclic property of traces and determinants, 
L{A) = L{P~^AP)^ and the latter evaluates to: 



d/Bou-t 


0 0 
0 ^/K 0 


Sa 

6a^ 

6N 


Cxdt 


"1 o' 

d,(3{n 

0 1 

dPtn 


DdBir. 


(5) 


L{A) = L{P-^AP) = 4|A + ti\‘^Re{p) (10) 

This will change sign from negative to positive when pass¬ 
ing through a Hopf bifurcation. 


Together, Eqs. ( [4p| ) may be written formally as: 


B. Scaling Laws 


dx = Ax dt + Bd/3[^ + Fdw (6) 

dPout = Cxdt P Ddpin (7) 

which is the standard form for a linear stochastic input- 
output system. 

Equation 0 separates the dynamics into three parts: 
a deterministic term Axdt, noise due to quantum fluc¬ 
tuations of the input B d^i^^ and additional free-carrier 
noise Fdw. (Here, dw is a vector Wiener process con¬ 
structed from the real and imaginary parts of the noise 
terms d/3^, djd^^ dPfca^ dwj, and normalized to satisfy the 
Ito table dwidwj = Sijdt; the matrix F is constructed so 
that Q is satisfied). 

The matrix A has three eigenvalues. Due to its 
doubled-up structure, complex eigenvalues must come in 
conjugate pairs. Thus, A can either have three real eigen¬ 
values or one real eigenvalue and one complex conjugate 
pair. If the equilibrium is stable, all three eigenvalues 
must have a negative real part. 

There are two ways for an equilibrium to go unstable. 
Eirst, a negative real eigenvalue can cross zero and turn 
positive. Since only a single direction goes unstable, the 
equilibrium point bifurcates into two stable equilibria. 
This is the standard cusp catastrophe of optical bistabil¬ 
ity in Kerr and cavity QED systems [26] . In our previous 
paper m, we discussed it in the context of carrier-based 
switches and amplifiers. By calculating the determinant 
of H, we can catch this instability - for stable equilib¬ 
rium, detH < 0, but if the equilibrium transitions to 
unstable, det A will become positive. 

Self-oscillation takes place when a conjugate pair of 
eigenvalues cross the imaginary axis. In this case, two 
directions go unstable, so the equilibrium point bifurcates 
into a ring of steady states, or more often, a limit cycle. 
The determinant will remain negative, but the product 


L{A) = (tr(H)^ — tr(H^)) tr(H) — 2det(H) (8) 

changes sign at this bifurcation. To see why, suppose 
that the matrix A has eigenvalues A, /i, /i*. Then for some 
transformation P, 


P-^AP = 


X 




(9) 


Equations ([l][^, and the resulting matrix H, have 8 
free parameters. That’s a lot. Naively, searching for os¬ 
cillating conditions would appear difficult because of all 
the parameters one must consider. However, several scal¬ 
ing laws let us reduce this to 6 “normalized” parameters, 
of which 3 are material constants. 

Start with equations of motion Ql- Let k = nPr] he 
the total cavity linear loss. Scale time, the electric field, 
the input field, and the carrier number as follows: 


t 


(y, —y 




Ai 


Ai 




N ^ 


N 


Intuitively, time t is scaled so that the cavity photon 
lifetime is one. The carrier number is scaled so that 
N = 1 shifts the cavity by one linewidth. The intra¬ 
cavity field a and input field /din are scaled to the two- 
photon absorption: |q:| = 1 means that the single- and 
two-photon loss processes are equally strong. 


Name 

Description 

CaAs PhC“ 

Si At-ring, 

k 

nPf] 

0.42 ps“^ 

0.31 ns“^ 

K 

I/O Coupling 

kl2'= 

k 

r] 

LA 

kl2 

0 

d 

TPA 

7.9 X 10“®fc 

3.7 X IQ-^k 

X 

Kerr 

Qd 

0 

(5 

FCD 

2.7 X IQ-^k 

(5.6-0.4i) X 10“^A; 

7 

Carrier Decay 

1.2k 

1.0k 

6 

62/5i 

0 

0.07 

c 

Slip 

34 

150 

X 

x/k 

0 

0 

7 

ilk 

1.2 

1.2 

R 

njk 

0.5 

1.0 

V 

r)lk 

0.5 

0 

A 

A/k 

varies 

varies 


“GaAs, huj = 0.9Eg, V = 0.25, Q = 5000, r/c = 2 ps; see fTol . 
compare m- 

^Si, A = 1.5)am, V = 40, Q = 4 x 10^, Tfc = 3 ns; see |19| . 

^All dimensional quantities in this table are scaled to the linear 
loss k. 

^Negligible, as dispersive effect is dominated by free carriers. 

TABLE I: Free-carrier cavity parameters for a state-of-the art 
GaAs photonic crystal and Si microring cavity. 
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Input 


Input 


FIG. 1: Oscillation region as a function of cavity parameters. 
Two materials are shown: Si at 1.5 fim. (left) and GaAs near 
the band edge (right). Oscillations occur to the right of the 
solid curves. Gurves represent different values of 7 , from 0.33 
to 3.0. Optical bistability occurs in the dashed region. Golor 
represents the steady-state input power. 


The reduced equations take the following form: 


da = [—[ 1/2 + 6N) — a'^a — i (^A ^ N)] a dt 

— \/R^indt + n 2) 

k 

dN = [7)^(a:*(a) + ~ 7 ^] dt + ^^^dw{12) 

In the absence of noise, these equations have 6 indepen¬ 
dent parameters: 


(5 


^1’ 






} 

} 

} 


Material 

Properties 

Cavity 

Design 

Tunable 


(13) 


where k = k, rj and Sc = Si — iS 2 . 

Once a material and laser wavelength are picked, only 
three parameters can be varied. The relative linear ab¬ 
sorption f] = \ — R typically cannot vary much - in a 
linear-absorption cavity it should be 0 ( 1 ) to maximize 
the nonlinearity, and in TPA materials like silicon it is 
zero. The ratio of optical to free-carrier lifetimes, 7 , 
can vary by several orders of magnitude, depending on 
the cavity geometry and Q. For instance, it is easy to 
make low-Q cavities with a very small 7 . State-of-the-art 
micro-rings have Q r\j 10 ^ and Tc ^ns and consequently 
7/fc r\j 1 . Coincidentally, photonic crystals tend to have a 
similar ratio, though the carrier decay mechanism (diffu¬ 
sion) is different. It is also possible to make large cavities 
with very high Q and large 7 . 


FIG. 2 : Oscillation region as a function of cavity parame¬ 
ters. Here, the x-axis is normalized input field rather than 
normalized N. 


Obviously, both the input power and detuning can also 
be varied. For a given material, these quantities exhaust 
the parameter space. By plotting the self-oscillating re¬ 
gions as a function of A and N (a function of the input), 
for reasonable values of 7 , we are essentially plotting the 
entire parameter space. As shown in Figure in a large 
fraction of the parameter space, the cavity should self- 
oscillate. 

Figure shows the self-pulsing region as a function of 
input field and detuning. This is generally similar to Fig¬ 
ure although the low -7 regions appear more accessible 
because, although the internal carrier number is high, the 
carriers are long-lived and the cavity requires less optical 
power. However, these cavities are complicated by opti¬ 
cal bistability (which occurs in the same region), and the 
slow response time is generally not desirable. The most 
desirable conditions seem to occur when the photon and 
carrier lifetimes are comparable, and the cavity is driven 
with a slightly detuned pump. 


II. SEMICLASSICAL SIMULATIONS 


Quantum simulations (in the semiclassical Wigner pic¬ 
ture) add noise to this model. For concreteness, in this 
section and the sections that follow, we consider a GaAs 
photonic-crystal cavity with parameters given in Table [l| 
however, our results are applicable to a range of devices. 
Quantities with units of time or inverse time (t. A, etc.) 
will be normalized to the cavity lifetime 1 /k. 

Figure shows simulations for a detuning A = —0.8. 
The input field is stepped from Pin = 25 (blue) to 175 
(black) in increments of 25. The top plot shows a typical 
time trace. Oscillations clearly set in at around Pin = 75. 
In addition to the amplitude, the oscillation frequency 
also increases with pump power. 

The right panel of Figure plots internal photon num- 
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FIG. 3: Top: time trace of Re[(a(t)] as the input field is 
stepped from /3in = 25, 50, 75,100,125,150,175. Bottom: 
Output field quadratures at these input powers. Right: Os¬ 
cillation between photons and carriers. 


FIG. 5: Left: Photon conversion efficiency, the ratio of limit 
cycle photons emitted to photons absorbed. Right: Ampli¬ 
tude gain /3out,a;/An,a; at cj = 1.7 for different values of seed 
amplitude /3in,a; = 2, 5,10. 



- 0.4 - 0.2 0.0 


0 3 6 9 12 15 18 21 24 27 30 


pump powers and detunings for this system. If a limit 
cycle forms, we are interested in its amplitude and fre¬ 
quency. The amplitude should be large, so that a sig¬ 
nificant fraction of the pump is converted to photons at 
the limit-cycle frequency. The frequency should be large 
enough that the pump and limit cycle fields can be eas¬ 
ily demultiplexed with a cavity. Figure plots both of 
these figures of merit. As expected, the amplitude at 
uj only becomes nonzero in the unstable region where 
^^[Amax] > 0. The frequency also grows with pump 
power, starting at uj ~ 1.7 and growing to u; ~ 4; this is 
probably a nonlinear effect of the strong pumping. 

Two other figures of merit are the limit cycle “effi¬ 
ciency” and the gain. Efficiency is defined in terms of 
the output and absorbed power: 


FIG. 4: Left: Stability of equilibrium point, measured by the 
real part of the largest eigenvalue of A. Right: Amplitude 
of limit cycle, with contours designating the limit cycle fre¬ 
quency. 


ber (horizontal) against carrier number (vertical). This 
provides a qualitative picture of the oscillations: when 
the photon number is high, more photons are absorbed 
and the free carrier number increases. Eventually the 
carrier number becomes so high that the cavity shifts off- 
resonance, reducing the cavity’s effective driving strength 
and consequently the photon number. Once the photon 
number falls, the carrier number falls because fewer pho¬ 
tons are being absorbed, but eventually this brings the 
cavity back on resonance, increasing the photon number 
and repeating the cycle. 

To get a more general picture, consider all possible 


^ _ -PuJjOUt 

-Puj,out -Pahs 

Efficiency is defined this way rather than output over 
input because much of the input power is not consumed 
by the device; it is just a constant bias that can be re¬ 
cycled. If there finite conversion to cj and no absorption, 
we say the efficiency is 1; if no conversion, it is obvi¬ 
ously zero. The left panel of Eigurej^ plots efficiency as 
a function of detuning and input field. While not close 
to 100%, the efficiency is not too small, either - peaking 
at around 20%. 

If we drive the device with a sinusoidal field whose 
frequency is close to the limit-cycle frequency, that field 
should be amplified. In this way, the free-carrier cavity 
acts as a phase-insensitive amplifier. The amplitude gain 
G(cj) = /3uj,out//^uj,in is plotted at CJ = 1.7 in the right 
panel of Eigurej^ 
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III. BELOW THRESHOLD: LINEAR 
AMPLIFICATION 


Below the Hopf bifurcation, a complex pair of eigen¬ 
values approach the imaginary axis. The corresponding 
eigenvectors span a plane in phase space; since motion 
tangent to this plane is only marginally stable, pertur¬ 
bations will be strongly amplified. Since this plane is 
two-dimensional, we expect linear, phase-insensitive am¬ 
plification of both quadratures of the input field HE!]. 

For any mesoscopic linear amplifier, an important 
question to ask is: how much noise does the amplifier 
have? Quantum mechanics sets a strict bound on the 
noise of a quantum linear amplifier m, and this bound 
is realized with the non-degenerate OPO [7]. Since free 
carriers are excited incoherently, one expects an amplifier 
driven by carriers to be noisier than a quantum-limited 
amplifier; however, if the difference is not too large, the 
free-carrier amplifier may still be preferred because of 
material, power, or footprint considerations. 


A. Nondegenerate OPO 


Although this paper is about free-carrier effects, we in¬ 
troduce the non-degenerate OPO here as a “benchmark” 
system because it is a well-studied system that saturates 
the Caves bound. It can be modeled as a quantum input- 
output system EUHH] with three fields: signal a, idler b 
and pump c. The internal Hamiltonian is 


H = A^a^a -f + AcC^c -h 


e'ahc^ — ea^h'^c 
2 i 


(14) 


and input-output couplings 


Li = ^/l^a 

L 2 = 

Ls = ^cC (15) 


Following the Wigner method of m , one can convert 
the master equation into a PDF for the Wigner function, 
and truncating higher-order terms, this PDE becomes 
a Fokker-Planck equation. This can then be converted 
into an SDE, and solving the SDE produces trajectories 
that sample from the Wigner function [ 22 ] . Adiabatically 
eliminating the pump field and setting A^ = —A 5 = A, 
Ka = K,b = (symmetric doubly-resonant cavity), one 
obtains the following equations of motion: 


dai 


da 2 


, 1 

d^out , 2 


—iA — 


■ [3 a2a2\ 

) 


“h e 0^2 


dt 


- \/^a2d/3in,3 


iA — 




0^2 + e 


dt 


-VKdl3in,2 - 
\fdaidt + d/?in,i 
y/na2dt -1- d/?in,2 


(16) 


(17) 

(18) 
(19) 


where = e*e//^ is the intrinsic coupling strength of the 
OPO. 

Here, ai and 0^2 are the signal and idler, which 
have the same lifetime but opposite detunings. The 
pump does not resonate. These equations are symmet¬ 
ric with respect to ai Because of the symme¬ 

try, the dynamics can be decomposed into a “symmetric” 
mode = {ai -h ol 2 )I‘^ “antisymmetric” mode 

a- = {ai — 01 ^ 2)12 (and likewise for the dp±). In addi¬ 
tion, define dwi^dw 2 as quadratures of the pump noise, 
dPin ,3 = {dwi + idw2)/2. The equations of motion be¬ 
come: 


da± = 


(—iA — k,/2 ± e)a± — ^ (a^ — 


2 \ 

a^)a^ 


dt 


-^/ddPin,± T {a±dwi - iazfdw 2 ) ( 20 ) 


The symmetric mode has gain (a +e term) while 
the antisymmetric mode has additional loss. As a result, 
at near- or above-threshold pumping, can become 
very large, but a- always stays near zero. In the weakly 
coupled case (^d <C 1 ), we can throw away the terms that 
couple and a- in the equation above, and combine 
the noise terms, giving: 


r P 2 ' 

daj^ = (—iA —/i :/2 + e)(a+— — |(a+| 


dt 


( 21 ) 


Linearizing about the fixed point q^i = 0^2 = 0, and 
transforming into the frequency domain, we arrive at the 
input-output relation: 


' (_(„_A) + iK/2f + (e/2)'i ■' 


+ 


cosh rj 

2 (K/ 2 )(e/ 2 ) 


(-(w - A) + iK/2) + (e/2)2 


( 22 ) 


sinh r] 


For phase-insensitive amplification, the gain G an noise 
S at frequency uj may be defined as: 
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0.5 1.0 1.5 2.0 2.5 0.5 1.0 1.5 2.0 2.5 0.5 1.0 1.5 2.0 2.5 


FIG. 6: Plots of the amplitude gain (top) and noise (bottom) 
for free-carrier cavity with A = —1.0 approaching the Hopf 
bifurcation. In the lower graph, the blue line comes from 
numerical simulation, the red curve is the analytic linearized 
model, and the black dashed curve is the Caves bound. 


G{u) 

S{co) 


^out,l {^) 
An.l(w) 

V2P(w), P{oj) 


(23) 

(5(w - u') ^ ^ 


In terms of 77 , they are: 


G{uj) = cosh 77 , 


S{u) = y 


2 cosh^ 77 — 1 


(25) 


with state x{uj) = , a*{—uj), N(uj)) and input- 

output field ;5(co’) = (;5(cc;),/3*(—cj)). This is the standard 
frequency-domain form for doubled-up variables m- 
Solving for x, this becomes a linear input-output rela¬ 
tion with a transfer function and a noise matrix: 


^out(^) 


D^C— 


—iuj — A 


rB 


An(w) 


+ 


r(a;)A„(c^) 

C ' 


—iuj — A 


w{uj) 


N {uj)w{uj) 


Applying the definitions of G and S in 
we find: 


Eqs. 


+ N{uj)N{lo)'‘ 
(29) 

Unlike the OPO, the free-carrier amplifier does not 
have a simple expression for G{uj) or S{uj). However, 
they are straightforward to evaluate numerically, and can 
be compared to a full nonlinear simulation. 

Figure !^ sh ows the gain and noise for the cavity studied 
in Sectiori|n| with A = —1.0. Far from the limit-cycle 
frequency, there is no gain and the output noise matches 
that of the vacuum. As the power is increased and the 
system approaches the Hopf bifurcation, the gain and 
noise at the resonance obviously diverge. But the noise 
always remains a factor of ^2 — 3 above the Caves 
bound (in terms of noise power, a factor of ^ 5 above 
the bound). This is due to the incoherent nature of the 
free-carrier nonlinearity. 


Gicu) = |T(w)n|, 


Sicof 


T{aj)T{uy 


Note that this S{uj) is different from the squeezing 
spectrum of l22l[2n]; rather, it is a measure of the elec¬ 
tromagnetic energy at frequency 00 . The squeezing spec¬ 
trum, by contrast, is a power spectrum of a homodyne 
measurement. 

From (25) one sees that the non-degenerate OPO sat¬ 
urates the Caves bound for phase-insensitive amplifiers 

HZ]: 


5(w) > y2G'(w)2 - 1 (26) 


B. Pree-Carrier Amplifier 

Turning to the free-carrier amplifier, first transform 
Equations ( [6][7| ) to the frequency domain: 


IV. NEAR THRESHOLD: CRITICAL 
EXPONENTS 

Near the bifurcation point, the system transitions from 
a stable fixed point to a stable limit cycle. Dynamical 
systems exhibit universal behavior near this bifurcation, 
in the sense that every system with a Hopf bifurcation 
can be transformed into the same normal form Ellsl]. 
The same is not true when one adds noise and quantum 
effects. Two systems with the same semiclassical equa¬ 
tions of motion can behave very differently once quantum 
noise is added. Nevertheless, all systems will show the 
same qualitative behavior near a bifurcation point. 

Before discussing the free-carrier oscillations, consider 
the non-degenerate OPO near threshold. Below thresh¬ 
old, there is a stable fixed point at a-y = a- = 0. Above 
threshold, there is a limit cycle at: 


— iujx{uj) = Ax{uj) ^ Fw{uj) (27) 

,5out(w) = C'S(w) + r>/3in(w) (28) 


l«+l = 


2e — tv 


(30) 










































FIG. 7: Left: Free-carrier limit cycles just above the bifur¬ 
cation point (noiseless simulation), for evenly spaced /3in = 
78,79,80,... Right: Size of the limit cycle in terms of a 
(blue) and |/3out(^<<^)|^ (black), and the critical exponents a ~ 
^out(^) ~ V^An 


Thus, if we smoothly vary the parameter e near the 
bifurcation point, e = /^/2 + (5e, the limit cycle amplitude 
goes as ^/e. This is a universal feature. However, not all 
OPO’s are equal up to a transformation - the behavior 
of the quantum states depends strongly on the value of 
p. For /3 ^ 1, dissipation is dominant and the system 
stays in a classical state with a positive Wigner function. 
For /3 ^ the Wigner formalism breaks down. (This is 
true for OPO’s in general. It is known that in this regime 
the degenerate OPO can access “highly quantum” states 
with non-positive Wigner function such as number states 
and cat states [3M34] .) 

The fixed-point eigenvalues near the bifurcation are: 
A = (e — /^/2) ± iA, and therefore: 


i«+i 


/Re[A] 

~W 




Re [A] 


(31) 


In classical dynamical systems theory, we can freely 
transform the system variable a, so the parameter [3 can 
be rescaled to 1. This is part of the process of trans¬ 
forming to the normal coordinate frame. Classically, a is 
dimensional and therefore [3 is not universal in any way. 
But in quantum mechanics, there is a universal scale for 
a\ the single-photon scale. Because of this, /3 becomes a 
universal parameter, and is related to the “quantumness” 
of the bifurcation. 

Figure shows that the free-carrier Hopf bifurcation 
satisfies the same critical exponent as the non-degenerate 
OPO: in terms of the input power An, the average oscil¬ 
lating field goes as \a\ ^ calculate the 


effective [3 for this bifurcation using Eq. (31): fitting to 
the figures, it works out to /3 ^ 0.0002, well in the semi- 
classical regime. 


FIG. 8: Phase plots (axes are Re[o], Im[a]) of the limit cycles 
for free carriers (A = —1.0,ain = 72.5 through 84.5) and the 
non-degenerate OPO (/3 = 0.0002, e = 0.48 through 0.52). 


Even after accounting for the free-carrier and OPO 
Hopf bifurcations are not equivalent up to a transfor¬ 
mation, as they would be in classical bifurcation theory. 
Again, the culprit is quantum mechanics: the incoherent 
process of carrier excitation and decay adds extra quan¬ 
tum noise, making the free-carrier limit cycle “fuzzier” 
than its OPO counterpart. This is shown in Eigurej^ 


V. ABOVE THRESHOLD: LIMIT CYCLE 


Above threshold, we classically expect a limit cycle. 
Quantum noise will blur this out to some degree, but 
sufficiently far above threshold, the cycle should be clear. 

Limit cycles are a classic topic in dynamical systems; 
some key results are reviewed in Appendix [A} To sum¬ 
marize the important points: Eor an n-dimensional phase 
space, there is a function that maps the 

limit cycle phase ^ and local perturbations u onto a por¬ 
tion of the phase space. When the perturbations are 
small compared to the limit cycle, they can be ignored 
entirely, reducing the dimensionality of the system from 
n to 1. This reduced system has the following equation 
of motion: 


d^ = Lodt + Y^ Re[Ri+ F{i)dw (32) 

i 


Here, is the response to an external perturbation 

dPin,i and F{^)dw is the intrinsic limit cycle noise. 

Any limit-cycle system can be used as a homodyne 
detector. To see why, consider a coherent input An,i = 

(A) where ujc is the limit cycle frequency. 

Averaging over many cycles, this input changes the limit- 
cycle phase as follows: 
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A^-cot = / VRe[Bi(e)*dAn,i] 

Jo i 

That is, the phase change has a normal distribution, 
with mean and variance given by the drift and diffusion 


constants: 

(R,(C)*e-*«)^ 

(34) 

Di = 


(35) 


(where ■ ■)dfj 



The drift term governs the response rate of the 
limit cycle to an external stimulus (in this case, the field). 
The diffusion term tells us how quickly the limit- 
cycle phase diffuses in the absence of a stimulus (assum¬ 
ing coherent inputs). Both terms show up in the homo¬ 
dyne measurement ( [^ . The standard quantum limit [7] 
bounds the accuracy of this measurement: in terms of the 
and this gives rise to a drift-diffusion inequality: 


+ [ F{f)dw 
Jo 

D^t] (33) 




(36) 


This relation holds for all limit cycles. One can also de¬ 
rive it from Eqs. ( [Mp5 ) by applying the Schwarz inequal¬ 
ity. Equality holds only for special, “quantum-limited” 
limit cycles where F{f) = 0 and Bi{f) ^ In the 

sections below, we compare the performance of the non¬ 
degenerate OPO and the free-carrier limit cycle using 
this metric, and show that the OPO saturates the drift- 
diffusion inequality, while the free-carrier device does not. 


A. Non-degenerate OPO 

Again, it will be important to contrast the results ob¬ 
tained here with the non-degenerate OPO; as we will 
show, this device can function as a quantum-limited ho¬ 
modyne detector for signal and idler fields. Because it 
is quantum-limited, no other limit-cycle device will beat 
the OPO at this task, just like no other linear amplifier 
can beat the non-degenerate OPO below threshold. 

As we show in Appendix [A} t he non-dege nerate OPO 
has a limit cycle with |<a+| = y^{2e — k)//3 and a phase 
that evolves as: 


df = A dt + Re 


= A dt + Re 


-iy/^ 

-iy/^ 

2a4 


d^i 


'in,+ 


dfd] 


'in,l 


2^1 ^ ’ 


(37) 



FIG. 9: Wigner function of the nondegenerate OPO {y = 
1.0,/3 = 0.01) subject to a bias /3i = O.lSz (red) and —O.lSz 
(blue). The state ^(t) for t > 0, which can be accurately 
read out with either homodyne or heterodyne detection, ef¬ 
fectively encodes a measurement of the p-quadrature of the 
input, Im[^i]. 


so that for signal and idler fields varying as /3ie 
the drift-diffusion terms are: 


= 

. Vk 

(38) 

*2K| 

Me,2 = 

*2K| 

(39) 

Di = 

K, 

(40) 

8|a+|2 


It is not difficult to see from ([^40) that the drift- 
diffusion inequality (36) is saturated. In this limit, the 


non-degenerate OPO functions as an optimal, quantum- 
limited homodyne detector. 

This is sketched in Eigurej^ Here, a non-degenerate 
OPO with A = 0 is used to measure the p quadrature 
of a signal field. Depending on the sign of the field, the 
state either drifts to the top or the bottom, and the dif¬ 
fusion incurred is due to the quantum uncertainty of the 
homodyne measurement. 


B. Pree-Carrier Cavity 

Since the equations of motion for the free-carrier cavity 
are more complicated, a simple analytic expression for 
and does not exist. However, these can be computed 
numerically. Eollowing the results of Section [Hlj it is rea¬ 
sonable to expect diffusion rates 5-10 times faster than 
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FIG. 10: Left: Limit-cycle phase diffusion for free-carrier cav¬ 
ity, A = —1.0, as a function of input field. Center: Phase dif¬ 
fusion for non-degenerate OPO, /3 = 0.0002, as a function of 
pump. Right: Combined, where the drift term 1 |/i^,i|^is 

the common x axis. 


for the non-degenerate OPO, the extra diffusion due to 
incoherent processes involving free carriers. 

Figure P!q| plots the simulated phase diffusion constant 
for both the OPO and the free-carrier limit cycle. 
As one approaches the bifurcation, the diffusion rate in¬ 
creases and diverges from the linearized result (35), solid 
curves in the figure. However, far from the bifurcation, 
the linearized model agrees with the full simulation for 
both the OPO and free carriers. 

To compare the OPO and free-carrier cavity on equal 
footing, the right panel of Figure IT plots the diffusion 
against the right-hand side of (36): ^ The 


OPO simulations, at least for large^a+j, lie on the line 
D^ = \ (green line), while the free-carrier sim¬ 

ulations lie a factor of ^ 10 above. 


C. Entrainment 


If the system is driven with a periodic seed field whose 
frequency does not exactly match the limit-cycle fre¬ 
quency cjc: Ihe limit cycle may or may not lock to the seed 
{entrainment)^ depending on its amplitude. To study this 
effect conceptually, assume a symmetric, noiseless limit- 
cycle model with a periodic drive and 

transform to comoving coordinates C = C “ oj[^t. Equa¬ 
tion (32) takes the form [6] : 


-^ = (cJc - CJin) - sin{C) (41) 

For frequencies \ujc — ^in| < |^An,a ;|5 there is a fixed 
point at C = sm~^{{ujc — ^in)/\Pin,ujB\), so the oscillator 
will lock to the seed. If we plot uji^ on the x-axis and 
on the y axis, this phase locking will happen in a 
vertical cone centered at (cJc^O)- Full free-carrier cavity 
simulations also show this effect. Figureshows results 



FIG. 11: Entrainment of free-carrier limit cycle, A = —1.0, 
An = 100. Top: Phase plots of the output field in a ro¬ 
tating wave frame, /3out (mean subtracted). For large 

seed inputs, the device clusters to one side of the diagram, 
indicating phase locking. Bottom left: output spectrum as a 
function of seed power, at cuin = 1.9. Bottom right: Entrain¬ 
ment cone. Plots of o(cjin) and a{(jUc) (intracavity amplitude 
at seed and natural frequency, respectively) as a function of 
seed frequency and amplitude. 


for a A = —1.0 cavity with pump An = 100, which natu¬ 
rally oscillates at ujc = 2.27. On top of this, an oscillating 
field drives the cavity. 

The top pane in Eigure shows the real and imagi¬ 
nary quadratures of the output field in a rotating-wave 
frame: . This is for seed frequency cjin = 1.9 and 

cavity frequency uJc = 2.3, so jcjin — uJc\ ~ 0-4, or about 
16%. Eor weak seed fields, the rotated output makes 
loops about the origin - the phase is not locked. How¬ 
ever, around An,a; = 10, A clusters in a given direction - 
indicating locking. 

The bottom-left plot shows the output spectrum 
Po\it{^) as a function of uj and the seed amplitude. One 
sees two peaks, one at the limit-cycle frequency uJc and 
one at the seed frequency cjin- The peak at the natural 
frequency uJc is strongest when the pump is weak, and 
eventually goes away for strong pumping. Conversely, 
the peak at the drive frequency cjin is absent for weak 
pumping, and grows with the pump strength. 

This is seen more clearly in the bottom-right plots. In¬ 
stead of confining ourselves to cjin = 1.9, in these plots 
we vary both the amplitude An,u; and frequency cjin of 
the pump. The left plot shows the power at the input 
frequency, while the right plot shows the power at the 
original frequency. Inside the entrainment eone, the os¬ 
cillator locks and the former dominates; outside the cone, 
the oscillator is unable to lock and the natural frequency 
is dominant. 
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FIG. 13: Optical free-carrier cavity used as a node in an Ising 
machine. 


FIG. 12: Left: Time traces of the limit-cycle phase ^ for a 
driven system where the seed phase jumps by one radian at 
t = 0. Right: Response rate l/r, obtained by exponential 
fitting, as a function of seed amplitude auj. Parameters: A = 

-1.0, Uin = 100. 


each spin maps onto an angle = (cosC^, sinC^) and the 
Hamiltonian becomes: 




(43) 


From the shape of the entrainment cone, we estimate 
B ^ 0.04 for this set of parameters. 


D. Impulse Response 

Suppose that the oscillator has been locked to an ex¬ 
ternal field and now the phase of that field is changed. 
The oscillator should follow that phase, but there will be 


a time lag. From Eq. (41) we can estimate this time lag 
to be of order: 


The general Ising problem for arbitrary Jij is NP-hard 

m- 

Ising problems map naturally onto oscillator networks. 
Let each Ising spin be mapped onto an oscillating free- 
carrier cavity. Let each oscillator have multiple indepen¬ 
dent input and output ports. This can be accomplished 
using the “railroad topology” of Figure Suppose that 
an output of cavity j is fed into an input of cavity i. As¬ 
suming all cavities have the same limit-cycle frequency, 
under the assumptions of Section [V^ the phase of cavity 
i evolves as: 




(42) 


dC,i = - Jij sin(Ci - Ci) 


(44) 


In Figure [T^ the same free-carrier system is simulated 
with a seed field cjin = uJc = 2.27. However, at time 
t = 0, the phase of the input shifts by 1 radian. For 
seed amplitudes ^ 3, the system quickly realigns 

to the new phase, with a time-constant given by (42). 
From this, we can estimate B « 0.02. This agrees with 
the entrainment-cone estimate to within a factor of 2; 
the lack of exact agreement is due to the circular cycle 
assumption that underlies (41, 42). 


VI. APPLICATIONS 


A. Ising Machine 


Many optimization problems can be recast as Ising 
problems, which involve finding the minimum of the Ising 
Hamiltonian: H = Jij^i • aj. If a is constrained to 
lie on the x^-axis the problem is called an XY model, the 


where Jij depends on the waveguide coupling, the phase 
of the connection, and the limit-cycle amplitude. It is not 
difficult to see that, with the appropriate connections, 
one can realize a cavity network that minimizes (43) by 
the steepest-descent method. 

A full discussion of optical Ising machines is beyond 
the scope of this paper. The concept was proposed by 
Utsunomiya et al. [1], who suggested implementing it us¬ 
ing injection-locked lasers. Recent theoretical work [36] 
and experiments with 4-bit m and 16-bit [38] Ising ma¬ 
chines using a time-multiplexed pulsed OPO show that 
the device matches or surpasses classical algorithms in 
accuracy. However, free-carrier oscillations may be a 
preferable platform for Ising machines because of their 
low power requirements and compatibility with existing 
fabrication processes. 

Figure shows the simulated Ising-machine perfor¬ 
mance for antiferromagnetic couplings on five graphs: 
pair, triangle, square, pentagon and tetrahedron. Of 
these, the pair and square have zero-energy configura¬ 
tions, while the rest are frustrated systems. The square 
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3-5 Spin Ising Solver 



achieve a per-watt performance orders or magnitude 
greater than a microprocessor solving the same problem. 
For the network used in Figure (see Sec. for cav¬ 
ity parameters), during oscillation each cavity consumes 
^ 2000 photons, or about 0.5 fJ, per cavity lifetime and 
takes ^ 100 lifetimes to converge, an energy cost of ^ 50 
fJ per spin and a computation time of ^ 300 ps. A micro¬ 
processor using steepest-descent or stimulated annealing 
will also take ^ 100 steps to converge, but be required 
to compute (44) at each step. Since (44) involves com¬ 
puting a trigonometric function, it will take ^ 50 flops 
and ^ 100 clock cycles per step [39], or ^ 5000 flops per 
spin overall. Presently, the most energy-efficient super¬ 
computer is the L-CSC at GSI, Darmstadt, which runs 
at 3 GHz and requires 0.2 nJ per flop [40], giving a sim¬ 
ulation time of ^ 3 /is and energy cost of ^ 1 /iJ per 
spin. On the basis of this rough calculation, the free- 
carrier Ising machine should perform ^ 10^ x faster and 
consume ^ lO^x less energy. 


FIG. 14: Ising machine performance for small graphs. Top to 
bottom: pair, triangle, square, pentagon, and tetrahedron. 


16-Spin Ising Solver 



FIG. 15: Ising machine performance for 16-gon and frustrated 
16-gon with cross-couplings. 


B. Free-Carrier Relay 

In a previous sections, we showed that free-carrier cav¬ 
ities can undergo spontaneous self-oscillation if driven 
hard enough. Here we show that this can be used to 
construct a free-carrier “relay”. Such a device has many 
logic applications, including message passing algorithms 
for error correction m- A relay acts like a classical 
GNOT gate: if the digital inputs A,Be {—1,1}, then 
the relay maps these to: 

{A, B) {A, AB) (45) 

That is, output B is flipped if A = — 1. 

The relay is a circuit with two free-carrier cavities, ar¬ 
ranged as in Figure [^ The inputs A and B arrive on 
the same channel, but are offset in frequency. Data is 
encoded on the phase of the inputs (0 or tt), not the 
amplitude; thus, for a fixed field amplitude |A|, a 1 cor¬ 
responds to +|A|, while —1 corresponds to — |A|. 

The input is mixed with a pump field on a beamsplit¬ 
ter, so that the field entering cavity a± is: 


and tetrahedron were studied with an OPO Ising machine 
in [37] . 

Larger networks also show convergence in reasonable 
time. In Figure [^ we plot the performance of a 16- 
spin network, both with a nearest-neighbor interaction 
and with a cross-interaction (which shows frustration). 
These are the graphs studied in the OPO network of [38] . 
As long as it does not get trapped in local minima, the 
device converges to the minimum of U[(] in 50 — 100 
cavity lifetimes. 

Because the free-carrier Ising machine maps the op¬ 
timization directly onto the hardware dynamics, it can 


An,± 


A±E^ 


(46) 


A free-carrier cavity will self-oscillate if the input field 
is stronger than some threshold: \/3i^\ > /3th- Let: 


1^1 - \Ep\ < /3th < |A| + \Ep\. (47) 

If A = -hi, then the top resonator is above thresh¬ 
old and self-oscillates at cj, while the bottom resonator 
does not self-oscillate. For the B field at this frequency, 
this means that the top channel has more gain than 
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FIG. 16: Left: Layout of the free-carrier relay. Right: Relay 
behavior when control bit A is set to +1 (left) or —1 (right). 


the bottom channel. When these are interfered on a 
beamsplitter, the output at this frequency is ^(Ghigh — 
Since Ghigh > Glow, the phase of B does 

not change. 

On the other hand, if A = —1, the lower channel has 
higher gain. When recombined on the beamsplitter, the 
output is — I (Ghigh — - the phase of B does 

flip. This is shown in Figure Thus, the relay realizes 

the CNOT map {A,B) {A,AB). 

Figure demonstrates the relay operation. Two re¬ 
sults are plotted: a “base” case with the same cavity 
parameters used elsewhere in the paper (blue in figure) 
and a hypothetical “lOx NL” case where the nonlinearity 
(parameters /3) has been increased by a factor of ten. 
Both cavities have a detuning A = — 2.0. In order to 
control the phase of the beam at cj, the input A must 
be fairly large (^4 = ±65 was used here, scaled by 
for the lOx NL case). However, the input B at cj can be 
quite small; in the simulation taking a value of about 3. 
Since the output amplitude is around 7, this provides an 
XOR with enough gain for a fanout of 4-5. 

Both relays display the same overall behavior, but be¬ 
cause the cavity in the lOx NL relay has a stronger non¬ 
linearity, it operates at a lower photon number and thus 
the photon shot noise is more significant. This degrades 
the performance of the XOR gate. Ultimately, there is 
tradeoff between gate fidelity and energy consumption 
for free-carrier based systems. Since this tradeoff arises 
from quantum mechanics, it cannot be avoided by choos¬ 
ing different materials or cavity designs. The benefit of 
our SDE approach OH is that it reveais not oniy the 
ciassicai behavior of the reiay, but aiso this basic quan¬ 
tum iimit to its performance. 


VII. CONCLUSION 

Systems with a Hopf bifurcation can perform a wide 
range of usefui tasks with appiications in sensing and 


FIG. 17: Left: Plots of the real and imaginary parts of the 
rotating-frame output B^, as a function of the input A and 
B^. Right: Time trace of the relay output (top), where the 
inputs A and B^ are switched regularly (bottom). Both the 
base (blue) and lOx NL (red) scenarios are shown. Outputs 
are scaled by for the lOx NL case. 


photonic logic. In this paper, we have studied the su¬ 
percritical Hopf bifurcation in a semiconductor optical 
cavity where the dominant optical nonlinearity is due to 
free carrier dispersion. Following our previous paper, we 
simulated the dynamics of a the free-carrier cavity using 
Wigner SDE’s that capture both the semiclassical mo¬ 
tion and the quantum fluctuations in photon and carrier 
number. 

Below the bifurcation, the free-carrier optical cavity 
acts as a phase-insensitive amplifier. This device is the 
basis for heterodyne detection, where both quadratures 
of the field are simultaneously measured with an added 
noise penalty. The Caves bound places a lower limit on 
the noise, and this limit is satisfied in the non-degenerate 
OPO. By contrast, the free-carrier cavity has ^ 5x more 
noise in the output, an effect we attribute to the incoher¬ 
ent nature of carrier excitation and decay. 

Above the bifurcation, the device has a limit cycle. 
Quantum fluctuations cause the phase of this cycle to 
diffuse, and the diffusion rate can be computed by lin¬ 
earizing the SDE’s in a normal coordinate frame centered 
on the limit cycle. In this limit, one can use the device to 
store a continuous number in the range [0,27r), or alter¬ 
nately, to perform a homodyne measurement on signals 
at the limit-cycle frequency. Limits on the efficiency of 
homodyne measurement lead to a quantum lower bound 
on the limit-cycle diffusion rate. This bound is saturated 
by the non-degenerate OPO, while the diffusion rate of 
the free-carrier cavity is ^ 10 x larger. Again, this is due 
to the incoherent carrier excitation and decay processes. 

Limit-cycle systems are useful in logic and computing 
because they can be locked to external signals, and their 
outputs can in turn be used to lock other limit cycles. 
While an analysis such large-scale networks is beyond 
the scope of this paper, we have explored the basic phe¬ 
nomenon that underlies this: entrainment in an external 
field. Utilizing entrainment, we showed that the free- 
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carrier cavity can be used to construct a coherent Ising 
machine that finds the minimum of a preprogrammed 
cost function. With reasonable cavity parameters, such 
a coherent Ising machine could run ^ 10^ x faster with 
^ lO^x less energy than a comparable algorithm on a su¬ 
percomputer. In addition, we showed that entrainment 
can be used to construct a limit-cycle “relay” - an all- 
optical classical CNOT gate, which has applications in 
message-passing schemes. 

Although the free-carrier cavity is noisier and performs 
more poorly than quantum-limited systems like the non¬ 
degenerate OPO, it is much more convenient to build. 
Free-carrier optical cavities can be built from silicon or 
III-V materials, which have mature and scalable fab¬ 
rication processes. In addition, the per-photon effect 
is much stronger, enabling operation at lower powers. 
When it comes to building an actual device, these prac¬ 
tical concerns may prevail over the theoretical elegance 
of quantum-limited systems. 


where uo is the oscillation frequency. The map x : M ^ 
defines the attractor’s manifold, and is sufficient if we 
are only interested in how the system behaves without 
forcing. However, the map tells us nothing about forcing 
or deviations from the attractor. When noise and forcing 
are present, these perturbations become relevant, and we 
need more information about the system to handle them. 


1. Linearization About Attractor 

Consider a nonlinear system of differential equations 
of the most general form: 


dxi = [fi{x) + Fi{x, t)] dt + gij{x)dwj{t) (A2) 
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Here, x is the state of the system, fi{x) is its natural 
(unforced) derivative, gij is the noise coupling (to Wiener 
process dwj) and Fi{x^t) is the external forci ng. I n the 
absence of forcing, let’s suppose that Equation ( |A2| ) gives 
rise to a stable attractor x{ujt). This has natural period 
T = 27r/co’, so ;r(^+2n7r) = x{^) for integers n. Deviations 
from this cycle are given by: x{t) = x{ujt) -\-Sx{t). In the 
absence of noise or external forcing, the perturbations 
evolve as follows: 


Appendix A: Limit Cycles and /c-dimensional 
Attractors 

Many dynamical systems do not have a fixed point. 
Instead, they have a stable limit cycle, or more generally, 
a stable /c-dimensional attractor. The k = 1 case corre¬ 
sponds to a limit cycle. The cycle may be parameterized 
as follows: 

x{t) = x{ujt) (AI) 



FIG. 18: Diagram of a limit cycle in a normal coordinate 
frame (left) and in the actual phase space (right), along with 
the transverse (blue) and longitudinal ViX (red) vectors. 


d{5xi) = ^^5xi = Aij{x^)Sxj (A3) 

OXj 


where Aij{x) = dfi/dxj is the Jacobian of the dynamical 
system; see 0, and ^ is the attractor phase, with x^ = 

x(0. 

The key trick is to perform a coordinate transformation 
that separates the dx, and n-dimensional vector, into I 
longitudinal perturbation and n — 1 transverse perturba¬ 
tions. The longitudinal perturbation keeps the system on 
the limit cycle, and therefore does not decay. The trans¬ 
verse perturbations deviate from the limit cycle, and de¬ 
cay to zero as t ^ oc. We denote these by and as 
follows: 


n—1 

5x{t) = + E “.(*)«£’ (") 

i=0 


Here we have traded an n-dimensional state vector x{t) 
for n — I transverse variables Ui (t) and one longitudinal 
variable S^i. 


Applying (A4) to the equations of motion with noise 
and forcing, we obtain: 
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d{S^)v^ + duie^^^ = ^ 


Ai- \ dvf 

A{x^)v^-lo — 




+ 




Dte 


(i) 


Dtv^=0 

Ui^dt + {F{x^^t)dt + g{x^)dw) 

(A5) 


(implicit summation over i) 

The covariant derivative Dt of a ^-dependent vector is 
defined as 


Dm = 


(A6) 


This derivative accounts for both the equations of mo¬ 
tion and our parameterization near the limit cycle. It is 
similar to the covariant derivative in Riemannian geome¬ 
try [42] . Because the tangent vector always transforms 
into itself when propagated around the manifold, its co¬ 
variant derivative is zero. Likewise, because the trans¬ 
verse vectors always decay to zero, they cannot evolve 
into v^; thus has no component. 


In matrix form. Equation (A5) is 


1 

'dismi 


r 1 

5^{t) 


du{t) 


0 Aej 



dt 

-\-F{x^, t)dt + g{x^)dw (A7) 
This becomes a matrix ODE: 


m 

u{t) 


n-1 r 


+ 






0 DfC^ 


m 

u{t) 


dt 


-1 


d^{t) = uj F Bl{ 0 {F{x^,t)dt + g{x^)dw) (All) 

du{t) = AT{^)u{t)dt F Bt{C) {B{x^^ t)dt + g{x^)dw) 

(A12) 

This equation captures our intuition regarding limit 
cycles and attractors. External forces (F, g) can give rise 
to two kinds of perturbations: longitudinal (encoded in 
changes to and transverse {u). Because of our choice of 
coordinates, the perturbations evolve independently. The 
At matrix causes transverse perturbations to decay as 
t ^ oc, while longitudinal perturbations do not. Often, 
we are only interested in the longitudinal perturbations; 
in this case we can ignore the u{t) altog ether. 

Altogether, we can arrive at ( |A11|A12 ) for an arbitrary 
limit cycle by following these four steps: 

1. Get equations of motion dx = [f{x) F F(x, t)] dt F 
g{x)dvu 

2. Get limit cycle x{^) and the tangent vector v^ 

3. Eind a set of vectors at each point ^ that satisfy 

the following: 

(a) spans the whole vector space 

(b) Perturbations along the Sx ^ eventually 
go to zero as t ^ oc 

4. Gompute At^Bl^Bt in Eqs. ( A8| A9 ) 


2. Non-degenerate OPO 

Now we apply this to the non-degenerate OPO intro¬ 
duced in Section El ^ The equations of motion are re¬ 
produced below: 


da 4 


{—iA — k/2 ± e)a± — ^ (^a'^a±a± — dt 

-^/^dPin,± F \ {a±dvui - (^13) 


0 0 

'mi 


'Bl 

0 At 

u{t) 


Bt 


{F{x^A)dt F g{x^)dw) 

{F{x^, t)dt + g{x^)dw) 
(A9) 


The limit cycle occurs at: 
(A8) 


i«+i = 


-k/2 

~JF~ 


(A14) 


In the equations above, ^(t) = ujt has a fixed time- 
dependence. The dynamical variable d^{t) adds a per¬ 
turbation to this We can roll 5^ into turning ^ into 
a dynamieal variable^ so the state vector becomes: 


n—1 


x{t) = x{^{t)) + y] ujiFF 

j=0 


(AlO) 


Eollowing the procedure above, we first find a mapping 
from [0, 27r] to the limit cycle. This is easy: (a+(C) = 
(a_(^) = 0. Next, one needs the v^ and In 
terms of the basis ((a+,a_), a good choice is: 




(A15) 

One can check that these are linearly independent (in 
doubled-up space) and span the whole space. Plus, due 


—ia^ 

, eW = 


, e(2) = 

0 

, e(3) = 

0 

0 


0 


1 


i 


The matrix ODE becomes: 
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to the symmetry of the problem, it should be pretty clear 
that perturbations orthogonal to the limit cycle or 

perturbations to the a- mode always decay to 

zero. 

In this case we are not concerned about deviations from 
the limit cycle, so there is no need to calculate the At 
(which depends on covariant derivatives D^e^). All we 
need to find is Bl- At the end of the day we get the 
following equation of motion: 


If the inputs An,i, An,2 are vacuum noise, the noise 
term on the right becomes 




e = A + Re 



(A16) 
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